Isolation, functional evaluation, and fermentation process optimization of probiotic Bacillus coagulans

Bacillus coagulans is a probiotic agent widely used in various industries. In this study, we isolated a novel strain of B. coagulans, X26, from soil and characterized its properties. X26 exhibited superior enzyme, acid, and biomass yields when compared with other bacterial probiotics and an antibiotic. Moreover, X26 significantly improved the body weight of rats, highlighting its potential for industrial development as a supplement for animals. To optimize the fermentation process of this bacterium, we adopted the response surface design. When X26 was cultured in a medium with 16.5 g/L maltose, 25.00 g/L yeast extract, and 3.5 g/L K2HPO4, the optimal yield was predicted to be 5.1 × 109 CFU/mL. Consistent with the prediction, the yield of X26 in a 500-mL flask culture was (5.12 ± 0.01) × 109 CFU/mL, and in a 30-L fermenter was (5.11 ± 0.02) × 109 CFU/mL, accounting for a 9.9-fold higher field than that with a basal medium before optimization. We further optimized the fermentation process in the 30-L and a 10-T fermenter, generating yields of (7.8 ± 0.2) × 109 CFU/mL (spore rate: 96.54%) and (8.7 ± 0.1) × 109 CFU/mL (spore rate: 97.93%), respectively. These yields and spore rates were achieved at 45–55°C, the typical fermentation temperature of B. coagulans. Our findings indicate that B. coagulans X26 is a promising probiotic with considerable potential for cost-effective industrial fermentation.


Introduction
In recent years, probiotics have attracted extensive research interest for applications in several industries, including medicine, food, agriculture, animal husbandry, and fisheries.Bacillus coagulans is a facultative anaerobic bacterium that can adapt to low-oxygen intestinal conditions and ferment lactic acid to produce L-lactic acid, which lowers intestinal pH and inhibits the proliferation of harmful bacteria.As such, B. coagulans also serves as a probiotic.Indeed, Zhou et al. introduced the application of B. coagulans in animal husbandry [1], whereas Mu and Cong highlighted the importance of B. coagulans in medicine and its therapeutic effects for acute diarrhea, irritable bowel syndrome, antibiotic-associated diarrhea, constipation, and colitis [2].Konuray et al. further revealed the widespread application of B. coagulans in the food industry as a probiotic [3].Accordingly, B. coagulans is considered to be a safe and effective probiotic.
The growth temperature of B. coagulans ranges from 30-60˚C; however, the optimal growth temperature is 45-55˚C, producing a yield of 1.15 × 10 9 -3.2 × 10 12 CFU/mL and a spore number of 3.38 × 10 8 -4.57× 10 9 CFU/mL [14][15][16][17].The spore yield obtained by fermentation at the typical B. coagulans growth temperature is significantly lower than that at a fermentation temperature of 30-40˚C.Hence, several factors limit the practical application of B. coagulans, especially its characteristics and performance under different fermentation conditions.Therefore, further screening of high-yield strains and optimization based on the biological characteristics of B. coagulans is needed to promote the industrialized production of B. coagulans under 45-55˚C culture conditions.
In this study, we aim to selectively isolate and screen B. coagulans strains with high levels of stress resistance from soil and establish an effective fermentation process and industrial development route to ultimately promote its commercial application.The findings of this study, specifically related to cultivation conditions, provide a useful reference for research in the fields of probiotic development and application.

Selective isolation and screening of B. coagulans
Fourteen mixed acid soils were collected from orchards and vegetable fields in different regions (S1 Table ); 3.0 g of each was weighed, mixed in 100 mL of sterile saline, shaken in a water bath at 37˚C and 150 rpm, and extracted for 1 h.The extract was treated in a water bath at 80˚C for 10 min, and the 5% inoculum was treated in a pH 2.0 yeast extract peptone-dextrose (YPD) liquid medium containing 0.3% porcine bile salts and 0.1 mg/mL 7012 pepsin (Sigma-Aldrich, St. Louis, MO, USA) at 37˚C for 2 h.Subsequently, the extract was pressed, transferred to a YPD liquid medium for enrichment culture over 18 h, and then diluted and spread in a Man-Rogosa-Sharpe (MRS) medium plate (pH 5.5) for isolation and purification.Colonies were selected for identification after culturing at 55˚C for 48 h.
16S rDNA was used for preliminary identification during the screening process.No further identification of non-target strains was performed.B. coagulans was first identified using 16S rDNA, and B. coagulans X26 was sent to the China Type Culture Collection Center-an International Depository Authority (IDA) approved by the World Intellectual Property Organization (WIPO)-for preservation and identification based on its appearance, physiology, biochemistry, and molecular markers.The deposit number of B. coagulans X26 is CCTCC-JD-2023061.

Soil collection statement
Currently, no regulated protocols exist for the sampling of soil from general locations, such as farmlands, in China.The soil samples used in this study were collected from vegetable gardens or orchards, and approximately only 100 g of each soil sample was required.Therefore, only verbal consent was obtained from the landowner for the collection of samples.

Growth and metabolic performance
A single B. coagulans X26 colony was transferred to an MRS medium and cultivated at a constant temperature of 45˚C for 24-48 h.The fermentation mixture was collected for lactic acid detection by high-performance liquid chromatography (HPLC).Xylanase production was determined using xylan as the carbon source as previously described [18].α-Galactosidase production was determined as previously described [19].

Bacillus coagulans spore detection
The bacterial solution was treated in a water bath at 80˚C for 10 min as the spore treatment condition, and viable bacteria were counted on an MRS medium plate after incubation at 45˚C for 36 h.

Ethics statement
This study was approved by the Experimental Animal Management and Use Committee of the Hubei Provincial Center for Disease Control and Prevention (project approval number 202120134).
Euthanasia was performed by cervical dislocation after anesthesia via intraperitoneal injection of 1% barbital sodium.

Liquid fermentation process
Shake flask culture.The glycerol tube was streaked on an MRS medium slant and cultured at 45˚C for 48 h.Single colonies were transferred to an MRS liquid medium and cultured on a shaker at 45˚C and 180 rpm for 8-10 h to form a seed solution.Next, 5 mL of a seed solution was inoculated into a 500-mL shake flask (loading volume 100 mL) and cultured at 45˚C and 180 rpm for 24-48 h.30-L fermentation tank process.We inoculated 100 mL of a seed solution into a 30-L tank (loading capacity 15 L) under the following culture conditions: 45˚C, pH 7.0, tank pressure 0.035-0.05MPa, aeration ratio 1.0 m 3 /(m 3 �min), stirring speed 100 rpm, and culture period 24-48 h.Foaming was controlled by the flow of silicone oil 1520.After feeding with 2 L of glucose and 2 L of yeast extract that was sterilized at 118˚C for 30 min and administered at a feed rate of 260 mL/h, fermentation conditions were regulated at a ventilation ratio of 1.6 m 3 / (m 3 �min), stirring speed of 200-400 rpm, and a dissolved oxygen rate between 10% and 50% which was maintained by adjusting the air volume and speed.
10-T production tank fermentation process.A total of 300 mL of seed solution was placed in a 500-L seed tank (loading capacity of 300 L) at a culture temperature of 45˚C, pH 7.0, tank pressure 0.03-0.05MPa, aeration ratio 1.0 m 3 /(m 3 �min), stirring speed 50-80 rpm, and culture period 8-10 h, and then transferred to a 10-T fermentation tank (loading capacity 6 m 3 ).The tank prefeeding process was conducted at 45˚C, pH 7.0, tank pressure 0.03-0.05MPa, aeration ratio 1.0 m 3 /(m 3 �min), and stirring speed 50-100 rpm.After feeding with 800 L of glucose and 800 L of yeast extract at a feeding speed of 104 L/h, fermentation conditions were regulated at a stirring speed of 100-150 rpm, ventilation ratio of 1.6 m 3 /(m 3 �min), and a dissolved oxygen rate between 10% and 50% which was maintained by adjusting the air volume and speed.
Preparation of live bacteria.The fermentation broth was cooled with a solid content of approximately 4% to 25˚C, concentrated approximately five times with a disc centrifuge, treated with 15% starch and 5% trehalose as a carrier, mixed well, and spray-dried (inlet air temperature: 160-185˚C; outlet air temperature: 65-85˚C) while stirring to prepare powdered B. coagulans X26.

Test design
Single-factor experiment.The basal medium (control treatment) comprised glucose 10 g/L, soybean meal powder 15 g/L, KH 2 PO 4 5 g/L, MgSO 4 �7H 2 O 0.2 g/L, and MnSO 4 1.0 g/L, at pH 7.0-7.2and was sterilized at 115˚C for 30 min.To the basal medium, 10 g/L glucose, 15 g/ L fructose, 15 g/L sucrose, 15 g/L maltose, 10 g/L lactose, 15 g/L maltodextrin, 15 g/L cornstarch, and 15 g/L bran flour were added as carbon sources for screening and comparison; 20 g/L soybean meal powder, 30 g/L corn steep liquor, 20 g/L yeast extract, 20 g/L peptone, 30 g/L corn gluten meal, 20 g/L corn peptone, soybean 20 g/L peptone, and 20 g/L yeast autolyzed powder were added as nitrogen sources for screening and comparison; and 0.2 g/L magnesium sulfate, 1.0 g/L manganese sulfate, 0.2 g/L calcium chloride, 3.0 g/L calcium carbonate, 3.0 g/L ammonium sulfate, 3.0 g/L dipotassium hydrogen phosphate, 3.0 g/L potassium dihydrogen phosphate, 3.0 g/L sodium chloride, 3.0 g/L sodium acetate trihydrate, 1.0 g/L L-glutamic acid, and 0.5 g/L cysteamine hydrochloride were added as inorganic salts or trace elements for comparison.
Plackett-Burman (PB) design.In the single-factor experiment, nine variables with a marked impact on B. coagulans X26 yield were screened and included as key factors in the PB design, including three carbon sources (lactose, glucose, and maltose), corn steep liquor, corn gluten meal, yeast extract, three nitrogen sources, manganese sulfate, dipotassium hydrogen phosphate, glutamic acid, and three trace elements.The PB design 0 level was determined according to the results of the single-factor experiment, and the high and low levels were set 1.25-2 times higher and lower than the 0 level, respectively.
Steepest slope test.Based on the PB design results, we optimized the horizontal concentration of the three most important variables to design the steepest slope test.The search direction was determined by the linear coefficient of the equation.According to the previous single-factor experiment, the increased step size was set as maltose 1.5 g/L, yeast extract 5.00 g/ L, and K 2 HPO 4 0.5 g/L; other variables were at the center point level.

Center combination design (CCD).
According to the PB design and the results of the steepest slope test, the CCD design was carried out with 16.5 g/L maltose, 25.00 g/L yeast extract, and 3.5 g/L K 2 HPO 4 as the center points at five levels (-1.68179, -1, 0, +1, and +1.68179).

Statistical analysis
Design Expert 11 software (Stat-Ease, Minneapolis, MN, USA) was used in the experimental design and for statistical analysis.The fit of model equations was determined using the coefficient of determination (R 2 ), the statistical significance of differences between groups was determined using Fisher's F test, and the significance of the regression coefficient was determined using Student's t-tests.The statistical significance of the data was judged by the p-value (< 0.05).

Isolation and screening of B. coagulans
To isolate novel B. coagulans strains with probiotic potential, we inoculated soil sample extracts on media at different culture temperatures and stress conditions.After screening 14 soil samples, two B. coagulans strains were isolated in the quadruple stress resistance test group incubated at 55˚C.Strain B. coagulans X26 was isolated from the MRS5.5 medium, whereas strain B. coagulans X60 was obtained from the nutrient agar (NA) medium.At 37˚C, no B. coagulans strains were isolated from any test group; bacterial species isolated under these conditions were primarily B. cereus [21] ).
The 12 B. coagulans strains, isolated from different sources, exhibited unique characteristics.The optimal temperature range of the 12 strains was between 45-55˚C; that of 4 strains was 45˚C, and that of 8 strains was 50˚C.Moreover, the 12 B. coagulans strains showed no growth at pH 3.5, rather, 11 strains grew slowly at pH 4.5, and 7 strains grew optimally at pH 7.0; pH range 5.5-8.5 had the optimal growth condition.All were facultative anaerobic strains that produced lactic acid; six strains produced xylanase, and another six produced α-galactosidase (S3 Table ).Although B. coagulans X26 grew well at 45-55˚C, its growth rate was slightly faster at 50˚C.Meanwhile, B. coagulans X26 grew considerably slower with a long lag period at 37˚C and did not grow at 28˚C (S1A Fig) .The optimal pH for B. coagulans X26 growth was 7.0 (S1B Fig); however, the strain tolerated 0.1 mg/mL 7012 gastric acid or 0.3% pig bile salt at pH 2.0.B. coagulans X26 produced lactic acid at a rate of 0.3655 g/g and also produced xylanase and α-galactosidase.
B. coagulans X26 exhibited advantages over the other strains in enzyme production, acid production, and biomass (S3 and S4 Tables).Due to its enzyme-producing ability, B. coagulans X26 may aid in host digestion; its production of acid may also improve the content of shortchain fatty acids to regulate the intestinal tract, while its high level of biomass may reduce the cost of its application.

Functional evaluation of B. coagulans X26
To further elucidate its probiotic effect, rats were fed B. coagulans X26, and their body weights were recorded weekly and compared to those in the oxytetracycline (a widely used antibiotic), B. licheniformis, and E. faecium treatment groups.During the feeding period, none of the rats had abnormal health condition.After feeding for 1 week, the body weights in the oxytetracycline, B. coagulans, and B. licheniformis groups were comparable to those in the control group, whereas those in the E. faecium group were significantly lower, which may be related to the gastric acid tolerance of E. faecium.However, after 2 weeks of feeding, the B. coagulans X26 group gained significantly more weight than the rats fed with other bacteria or the antibiotic, thereby indicating that B. coagulans X26 has a significant positive effect on rat growth (Table 1).As previously reported, B. coagulans has characteristics of Bacillus and Lactobacillus species [22] in that it is highly resistant to high temperatures, gastric acid, and bile salts, and produces lactic acid, which may explain its superiority in promoting rat growth.
Based on its superior performance in terms of biomass, acid production, enzyme production, and probiotic function, B. coagulans X26 was selected to optimize the fermentation process, which is expected to reduce production costs and facilitate its broad application.
PB design key factor screening.The PB experimental design and results are shown in Table 2.We used a 11-factor analysis table to perform variance and code equation analysis.The analysis of variance is presented in Table 3.The model F value was 9.32, which represents significance.Model terms with p < 0.05 were significant; factors C, D, G, H, J, and K were identified as significant model terms.Eq (1) shows that the response value caused by factors G, C, and D had the largest fluctuation range, and the prompt has the greatest impact on the response value.Factors G, C, and D were set as the key analysis factors.
The response surface in Fig 1A was relatively steep, indicating that B and C had a significant effect on the number of viable bacteria, and its contour map was approximately elliptical, indicating that their interaction was also significant.Based on the contour and response surface maps, at a maltose concentration of 16.5 g/L, the number of viable bacteria first increased and then decreased with increasing yeast extract content and increased with decreasing K 2 HPO 4 content.
Moreover, the response surface in Fig 1B was relatively steep, indicating that A and C had a significant effect on the number of viable bacteria; however, the contour map was not elliptical, thereby indicating that their interaction was not significant.The contour and response surface maps show that at a yeast extract content of 25 g/L, the number of viable bacteria increased with increasing maltose content and decreasing K 2 HPO 4 content.
The response surface in Fig 1C was relatively steep, indicating that A and B had a significant effect on the number of viable bacteria, and its contour map was not elliptical, indicating that their interaction was not significant.The contour and response surface maps show that at a 3.5 g/L K 2 HPO 4 content, the number of viable bacteria cells increased with increasing maltose content and decreasing yeast extract content.The quadratic surfaces of A, B, and C had extreme values, and the openings of the response surfaces were all downward, indicating that the fermentation formula could be optimized through this response surface optimization test to obtain the highest effect of viable bacteria count.
We analyzed the variance in the regression model (Table 5) and found that (1) the regression model F value was 278.06 and p < 0.0001, representing extreme significance.(A) Effect of the interaction of yeast extract and K 2 HPO 4 on the number of viable bacteria when the amount of maltose is zero (16.5 g/L).(B) Effect of the interaction of maltose and K 2 HPO 4 on the number of viable bacteria when the amount of yeast extract is zero (25 g/L).(C) Effect of the interaction of maltose and yeast extract on the number of viable bacteria when the amount of K 2 HPO 4 is zero (3.5 g/L).https://doi.org/10.1371/journal.pone.0286944.g001value was 1.88%, representing a high confidence level for the model, which can be used to analyze the change in the response value.
For the p-value, the primary terms A, B, and C, the interaction term BC, and the quadratic term B 2 had significant effects on the number of viable bacteria (p < 0.01); the effects of the quadratic terms A 2 and C 2 were also significant (p < 0.05).Moreover, the F value of each test factor reflected its degree of influence on the response value; with an increasing F value, the effect of the test factor on the response value increased.Given that F A = 57.08,F B = 562.14, and F C = 1392.81,the order of influence for these three factors on the number of viable bacteria was K 2 HPO 4 (C) > yeast extract (B) > maltose (A).
The above analysis indicated that the model-predicted value maintained a high degree of consistency with the actual value, indicating that the model was successfully established and statistically significant.The model was analyzed using Design Expert 11, and the optimal ratio was obtained as maltose 18.0 g/L, yeast extract 20.00 g/L, and K 2 HPO 4 3.0 g/L; the highest predicted yield was 5.1 × 10 9 CFU/mL.

Liquid fermentation verification and process optimization
The yield under the optimized culture conditions was verified by 500-mL shake flask fermentation.The number of viable bacteria measured after 36 h of culture was (5.12 ± 0.01) × 10 9 CFU/mL, which was consistent with the prediction results of the model, confirming good model accuracy and validity.Compared with that in the basal medium, the number of viable bacteria in the optimized medium increased by 9.9 times.A 30-L experimental tank was used to verify the CCD-optimized medium, and the measured viable bacteria count was (5.11 ± 0.02) × 10 9 CFU/mL, which was consistent with the 500-mL shake flask results and model prediction.
In the early stage of fermentation, the cells grew rapidly, dissolved oxygen levels decreased rapidly, OD 600 levels increased, and pH first decreased before it rebounded (Fig 2).Bacterial autolysis occurred in the later stage of fermentation, and the spore rate was relatively low.To prevent bacterial autolysis, we added glucose and yeast extract and controlled the dissolved oxygen level at 10-50% in the late logarithmic phase of bacterial growth before the stable phase, which increased the number of viable bacteria to (7.8 ± 0.2) × 10 9 CFU/mL.The spore rate reached 96.54%, and the number of viable bacteria increased by 1.53-fold before the adjustment.According to the culture differences between the 10-T production tank and the 30-L fermenter, we adjusted the process conditions and guaranteed parameters and found that the number of viable bacteria reached (8.7 ± 0.1)×10 9 CFU/mL with a 97.93% spore rate, which was 1.12-fold higher than that in the 30-L experimental tank.The fermentation broth had a solid content of approximately 4% and was cooled to 25˚C, concentrated approximately five times using a disc centrifuge, mixed with 15% starch and 5% protective agent as a carrier, and spray-dried while stirring to prepare a powdered B. coagulans X26; wherein the number of viable bacteria was (1.40 ± 0.1) × 10 11 CFU/mL, and the spore rate was > 98%.

Discussion
In this study, we screened B. coagulans X26 from soil under varying culture conditions and found that it showed a high level of resistance to multiple stressors, strong enzyme and acid production capacity, and high biomass yield.Meanwhile, B. coagulans X26 also improved the body weights of rats and the growth performance of laying hens [23], thereby highlighting its potential for applications in commercial animal husbandry.
Response surface design is a widely used optimization method.In fact, Chen et al. applied this method to increase the yield of B. subtilis WHK from (2.7 ± 0.75)×10 9 to (1.52 ± 0.06)×10 10 CFU/mL [24].Additionally, Sun et al. adopted the response surface design to optimize the fermentation conditions of Pichia pastoris GS115 to increase the production of xylanase.After optimization, the production increased by 2.1-fold [25].In the present study, we utilized response surface design to determine that the main nutritional factors affecting the X26 yield were yeast extract, maltose, and K 2 HPO 4 .After optimization of the medium, the yield increased by 9.9-fold.After process adjustment, the yield further increased by 1.53-fold.These results indicated that the yield improvement through process adjustment was lower than that of the previous medium optimization.However, adjustment of the post-fermentation process had a positive effect on improving the spore rate.In this study, at a fermentation temperature of 45˚C, the viable count in the 500 mL shake flask was (5.12 ± 0.01) × 10 9 CFU/mL, while that in the 30-L fermenter was (7.8 ± 0.2) × 10 9 CFU/mL, the number of spores was 7.53 × 10 9 , and the spore rate was 96.54%.Similarly, within the 10-T fermenter, the number of viable bacteria was (8.7 ± 0.1) × 10 9 CFU/mL, the number of spores was 8.52 × 10 9 CFU/mL, and the spore rate was 97.93%.Meanwhile, Dong et al. fermented B. coagulans BLCC1-0517 at 42˚C for 48 h using an optimized medium, and achieved a spore density of 3.38 × 10 8 CFU/mL [26].Li et al. cultured B. coagulans CMCC337209 in optimized conditions at 45˚C for 96 h to obtain a maximum viable count of 1.15 × 10 9 CFU/mL and a spore count of 1.09 × 10 9 CFU/mL [14].Sun et al. used high-density fermentation to culture B. coagulans CGMCC NO.7431 and generated a viable bacteria number of 3.095 × 10 9 CFU/mL at 45˚C for 16 h, although spore yield was not reported [15].Ding used high-density fermentation to optimize the yield of B. coagulans CNJD and generate a viable bacteria number of 3.2 × 10 12 CFU/mL; however, only 4.3 × 10 9 CFU/mL spores were obtained, corresponding to a spore rate < 1% [16].Qiu et al. evaluated the fermentation characteristics of B. coagulans T-8 cultured at 47˚C for 28 h and obtained a viable count of 4.8 × 10 9 CFU/mL and spore count of 4.57 × 10 9 CFU/mL [17].Hence, at a fermentation temperature of 45˚C, the spore amounts obtained in the current study in the 30-L (7.53 × 10 9 CFU/mL) and 10-T fermenter tanks (8.52 × 10 9 CFU/mL) were the highest reported values for B. coagulans fermentation at this temperature (S5 Table ).
In industrial production, the efficient recovery of live bacteria at low costs is the main challenge when developing high-content live bacteria preparations; two processes are commonly utilized.The first is freeze-vacuum drying which can be used for live bacteria or spores; however, its costs are high, thereby limiting its application in the animal husbandry industry.The second is the spray drying process, which is simple to perform and cost-effective.However, the high temperatures during the drying process kill cells and result in a low bacterial yield.To address this problem, we optimized the fermentation process to achieve a spore rate > 95% and considered the temperature tolerance of the spores to avoid extreme cell losses during the drying process.The screening of appropriate carriers and protective agents can further improve the dry yield of viable bacteria to 90%.Indeed, we prepared (1.40 ± 0.1) × 10 11 CFU/ mL B. coagulans powder, thus, confirming that B. coagulans X26 not only has probiotic properties, but also an appropriate yield capacity for industrial fermentation.

Conclusions
Our study describes an optimized industrial development protocol for B. coagulans X26, including strain screening, performance evaluation, effect evaluation, and 10-T fermentation.Herein, we provided an empirical basis for the subsequent isolation of B. coagulans strains from the environment and optimization of culture conditions.Simultaneously, we have highlighted the importance of strain type as a factor that influences the fermentation temperature and biomass of B. coagulans.However, based on the limited number of B. coagulans strains and fermentation conditions, this study may also have certain limitations.Further studies are warranted to assess the relationship between growth temperature and fermentation yield.
(2) R 2 = 0.9960 and R 2 Adj = 0.9924 represent a strong correlation, good model fit, and small experimental error.(3) The lack of fit term (0.1308) had a p > 0.05, representing good model stability, a small proportion of abnormal error in actual fitting, and overall good model fit.(4) The CV

Fig 1 .
Fig 1. Contour plots and response surface plots.(A) Effect of the interaction of yeast extract and K 2 HPO 4 on the number of viable bacteria when the amount of maltose is zero (16.5 g/L).(B) Effect of the interaction of maltose and K 2 HPO 4 on the number of viable bacteria when the amount of yeast extract is zero (25 g/L).(C) Effect of the interaction of maltose and yeast extract on the number of viable bacteria when the amount of K 2 HPO 4 is zero (3.5 g/L).

Fig 2 .
Fig 2. 30-L fermentation process.(A) 30-L fermenter process without feeding.(B) Feed process.https://doi.org/10.1371/journal.pone.0286944.g002 , B. siamese, B. amyloliquefaciens, B. veles, B. subtilis, B. chiella, B. licheniformis, B. paralicheniformis, and B. desert Sonora.B. stuarteri was obtained from different stress resistance test groups at 60˚C.These results suggest a low overall abundance of B. coagulans and that multiple stress treatments (80˚C for 10 min, pH 2.0, 0.1 mg/mL pepsin, and 0.3% pig bile salt) combined with a growth temperature of 55˚C are required to enrich B. coagulans in the soil bacterial community.These results also indicate that B. coagulans is extremely resistant to various stresses (S2 Table